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Abstract 

Photon emission arising in the course of transition between the states 
of quark-gluon and hadron plasma has been considered. Single-photon dis- 
tributions and two-photon correlations in central rapidity region have been 
calculated for heavy ion collisions at high energies. It has been found that 
opposite side two-photon correlations can serve as a sign of sharp transitions 
between the states of strongly interaction matter. 
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1 Introduction 



Forty years ago E.S.Fradkin and his students had calculated the photon po- 
larization operator in relativistic plasma at finite temperatures Q. These 
results will be used here for estimation of a new specific mechanism of pho- 
ton production which may appear effective for identification of transitions 
between the states of quark and hadron matter in heavy ion collisions. 

The phenomenon under consideration is the photon production in the 
course of evolution of strongly interacting matter. Let us consider photons 
existing in the medium at initial moment to having momentum k and energy 
u) in . Let the properties of the medium (its dielectric penetrability) change 
within time interval 5t so that the final photon energy is Uf. As a result of the 
energy change the production of extra photons with momenta ±k takes place 
these photons having specific two-photon correlations. Analogous processes 
were considered for mesons [§, ||, £|, [| and applied to pion production in 
high-energy heavy ion collisions |J. 

The conditions for a strong effect are the following: first, the ratio of the 
energies u in /u)f must not be too close to unity and second, the transition 
should be fast enough. 

2 Basic formulation 

Time evolution of the transverse photon creation and annihilation operators 
is given by canonical Bogoliubov transformation |7[] which represents solution 
of the Hamilton equations and contains two modes with momenta ±k : 



(polarizations are omitted for a moment). Here Bogoliubov coefficients u,v 
satisfy equation 



preserving canonical commutation relations and the limit t — > oo must be 
taken. Physically the process under consideration is analogous to parametric 
excitation of quantum oscillators. It was considered in more details earlier 



a 



o(k, t) = u(k)a (k, 0) +v(k)a f (-k,0), 
• + (-k, t) = u(k)a(k, 0) + w(k)a f (-k, 0), 



(1) 




(2) 
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The Bogoliubov coefficients u(k), v(k) are taken to be real valued and k 
dependent. So we use a parametrization 

u(k) = coshr(fc), v (k) = sinhr(fc) 



(3) 



thus introducing evolution parameter r(k). 

To get feeling of the main features of the evolution effect (and for further 
references and comparison) let us formulate a simple model - fast simulta- 
neous break-up of large homogeneous system at rest ||. In this case the 
resulting single-particle momentum distribution can be written in a simple 
form 



dN 
~d?k 



(at(k)a(k))| 



V 



t— >oo 



(2tt) 



n(k) cosh(2r(fc)) + sinh 2 r(k) 



(4) 



(for single polarization) where V is the volume of the system and n(k) is the 
level occupation number at t — 0. The first term in rhs of Eq.(4) describes 
amplification of existed particles and the second term describes the contri- 
bution arising due to rearrangement of the ground state of the system in the 
course of the transition. 

The transition effect is better seen in particle correlations. Two-particle 
inclusive cross-section is given here by 



d 2 a 



a d 3 k\d 3 k2 



(a{a 2 aia 2 ) = (a\ax) (c4a 2 ) + (c4a 2 )(a 2 ai) + (aja^a-ic^) (5) 



The fist term in rhs of Eq. (5) is the product of single-particle distributions, 
the second term gives the usual Hanbury Brown-Twiss effect (HBT) and 
the third term is essential if time evolution takes place giving opposite side 
photon correlations (see below). The correlators in Eq.(5) in the case under 
consideration have the form: 

(at(kx)a(k 2 )) = \n{k) + (2n(k) + 1) sinh 2 r(k)} t^G^x - k 2 ), (6) 



(27T)* 



(o(ki)6(k 2 )) = sinh2r(fc) 



n(k) + - 



V 



(27T)" 



rC(k 1 + k 2 ) 



(7) 



where G(ki ± k 2 ) represents normalized Fourier transform of the source vol- 
ume at break-up stage (G(0) = 1). It is sharply peaked function of ki ± k 2 
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(at zero momentum) having characteristic scale of the order of inverse size 
of the source, this scale being much less than characteristic scales of photon 
momentum distribution n(k) and evolution parameter r(k). So the last two 
functions may be evaluated at any of momenta k l5 k 2 ~ ±k (we suggest that 
the process is k — > — k symmetric). 

Relative correlation function which is measured in experiment is now 
given by 

C(k 1; k 2 ) = 1 + G 2 (k! - k 2 ) + R 2 (k)G 2 (k! + k 2 ) (8) 

with 

sinhr(/c) coshr(k)(2n(k) + 1) 
( ' ~ sinh 2 r(k)(2n(k) + 1) + n{k) ( ' 

As it can be seen from Eqs.(8-9), HBT effect is given simply by the form- 
factor G(ki — k 2 ) in this model whereas the transition effect depends strongly 
on evolution parameter r(k). In turn r(k) depends on time duration 5t of 
the transition. For very small characteristic times 5t the expression for r(k) 
is universal 0, 

rw =HsS)" ^ r<<1 (10) 

where u>i n and Uf are particle energies before and after the transition. For 
larger 5t the evolution parameter lessens. In general we expect that it falls 
exponentially at large u5t if the time dependence of the energy in the course 
of transition has no singularities at real times. So for large uj5t we shall use an 
exponentially falling expression motivated by solvable model expression [|]. 
Below, after necessary modification, we apply the above consideration to 
photon production in heavy ion collisions. 



3 Photons in plasma 

Spectrum of photons in plasma is given by dispersion equation 

ujl = k 2 + U(uj k ,k,T,iJ,,m) (11) 

Here II is the polarization operator for transverse photons dependent on 
temperature T, chemical potential fi and the mass m of charged particles. 
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Below we use an approximate form extracted from original expression JT| : 

uj 2 - k 2 ( ' ( oo + vk\ 



with 



k 2 \ \u — vk / 

2 />no / ^-»2\ ^/^ 



(12) 



^2 = 4gaT / d x ( x i _™\ nF ( x p/T) (13) 

7T Jm/T \ 1 Z J 

where a = 1/137, v 2 is the averaged velocity squared of the charged particles 
in the plasma, factor g takes into account the number of the particle kinds 
and their electric charges (g = 5/3 for u, d quarks) and tif is the occupation 
number of the charged particles (Fermi distribution). Polarization operator 
for scalar charged particles is approximately a half of that for fermions with 
substitution of Bose distribution for Fermi distribution. Evidently the po- 
larization operator plays the role of (momentum dependent) photon termal 
mass squared m 2 . 

We calculated the polarization operator and photon spectrum for three 
possible kinds of plasma: quark-gluon plasma (QGP) with u, d light quarks, 
constituent quark (m = 350Me^)-pion plasma and hadronic (pions and 
nucleons) plasma. Chemical potential (baryonic one) was taken to be equal 
to lOOMeV per quark corresponding to typical value for SPS energies. The 
temperature was taken to be equal to 140MeV (see below). 

The evolution parameter r(k) for photons is determined through photon 
energy u(k). For small momenta k the parameter r(k) is well approximated 
by simple expression (to be used for k? < AQMeV) 

, m 2 -, — 777^9 5m 2 , 
r(k) = 2* = —1^ kSr < i 14 

where mL are photon termal masses squared at both sides of the transition 
and (m 2 ) is their average mass squared ,cf Eq.(lO). At k = the values of 
5m 2 are equal to 289, 178 and 106 (in MeV 2 units) for QGP-hadron, QGP- 
valon and valon-hadron transitions correspondingly. Corresponding values of 
zero momentum evolution parameter r(0) are 0.330, 0.154 and 0.178. 

Higher momentum behaviour of r(k) (to be used for > 40MeV) is 
taken in the form 

6 2 

r (k) = ^exp(~k5r S j (15) 
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where 5t is time duration of the transition. The Eq.(15) is a simple version 
of the expression given by the solvable model 0] which is sewed together with 
Eq.(14) giving a monotonically decreasing function of the momentum. Below 
Eqs.( 14-15) will be used for estimation of the transition effect in heavy ion 
collisions. Only QGP-hadron transition will be calculated. In view of fact 
that evolution parameter r(k) appeared to be small number at all momenta 
k, all expressios will be taken in the lowest order in r(k). 



4 Transition effect in heavy ion collisions 

Let us now apply the above consideratins to photon production in heavy 
ion collisions. Let us suggest that the quark-gluon plasma is formed at the 
initial stage of the ion collision. Let the plasma undergoes expansion and 
cooling. The expansion is taken to be longitudinal and boost invariant ||. 
Recent lattice calculations || show rather low critical temperature of the 
deconfmement and chiral phase transition, T c « 150MeV as well as sharp 
drop of the pressure when the temperature approachs T c thus provocing 
instability in the presence of overcooling. So we do not expect long-living 
mixed phase and consider fast transition from quark to hadron matter with 
characteristic transition proper time duration 5t of the order of lfm/c. 

To calculate the transition effect one must shift to rest frame of each 
moving element of the system and integrate over proper times and space- 
time rapidities of the elements. Then single-photon distribution in central 
rapidity region y — reads: 



dN 



J rdr J drj J d 2 x T (po~^] 



d 2 kTdy 
+ I drj 



T 4- 7"W 
Iqgp + hr 



2 2pr c 2/ \ 
d {P) 



(16) 



with p = kx cosh 77. 

The first term in rhs of Eq.(16) describes photon production from hot 
quark-gluon plasma. Here _R 7 is the QGP production rate per unit four- 
volume in the rest frame of the matter iTOl: 



Po 



dR~, 



5aa< 



d 3 p 187T 2 



T 2 exp(-p/T) ln(l + 



Kp^ 

T ' 



(17) 
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with a = 1/137, a s = 0.4, k — 0, 58. It can be used also for hadron gas 
as its uncertainity is larger than the difference between the first-order QGP 



and hadron gas production rates [ IT| . Contribution from hadronic resonances 
are not considered here. The second term in rhs of Eq.(16) describes photon 
production due to transition from QGP to hadrons in the vicinity of proper 
time r c , cf Eq.(4). The time duration of the transition is taken to be small 
in this term in comparison with total time duration of photon production 
process. 

The photon production rate in Eq.(16) can be expressed through photon 
occupation number n(k): 

dR^ 2kr dn{k) 
Po l^ = J2^— (18) 

(with two polarizations included). Taking Eq.(18) into account one can see 
that if the velocities of the volume elements, as well as proper time interval 
in the first term in rhs of Eq.(16) are small then Eq.(16) is reduced to Eq.(4) 
as it should be. The photon occupation number n(k) in Eq.(18) appears to 
be small numerically, 

n(k) < 1 

That means in particular that transition radiation is dominated not by the 
photon amplification but by the ground state rearrangement (cf Eqs.(4-8)). 

As the last step one must specify temperature evolution. We suggest 
that the temperature depends on proper time of the volume element with 
power-like dependence: 

(T/T ) = (r/r )- 1 / b (19) 

where r and T are initial proper time and initial temperature. For final 
estimation we use 6 = 3 typical for hydrodynamical picture and choose low 
transition temperature T c = 140MeV. After transition the photons live some 
time in hadronic medium and we suggest termal momentum distribution of 
the hadrons (modified by the expansion of the system). We neglect ter- 
mal photon production below T c and do not introduce a special freese-out 
temperature. 

Below we will be interested in rather low photon transverse momenta 
kx (up to 500MeV) where transition effect is more pronounced. In this 
momentum region the QGP production term Iqgp depends mainly on final 
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temperature T c . For two main variants of initial conditions used in liter- 
ature [0: t Tq = l,T /T c = 3/2, t = lfm/c and r T = 1/3, T /T c = 
5/2, r = 0.2 fm/c a variable factor in /qgp changes inessentially (from 5.06 
to 4.37 for 6 = 3) and we will use for this factor an average value 4.7 in 
our estimations. The transition proper time r c also changes inessentially 
for these two variants of initial conditions being 3.00 fm/c and 3.02 fm/c 
corresp ondingly. 

The transition contribution 1^ in Eq.(16) appears essential only at very 
small momenta kr- So dealing with single-photon distributions one can use 
Eq.(14) for evolution parameter r(k). The resulting relative strength of tran- 
sition radiation 

R 1 {k T )=I§ ) /I QGP (20) 

appears sizable only in the momentum region kx < (15 — 20)Mev indepen- 
dently of the time duration of the transition. Ri reachs 4.44 at kx = and 
falls down to 0.06 at kx = 40MeV. Because of high background effects in 
this momentum region the single- photon transition effect should be difficult 
to observe experimentally. 

Mach better the transition effect is seen in photon correlations (cf Eqs. (5- 
9)) where it is first order effect with respect to r(k). Let us note that HBT 
effect for photons has now more complicated form than that in Eq.(8) because 
of finite time duration |13|] of the process of photon emission from QGP and 
it will not be exposed here. We consider only the transition effect (the 
third term in Eqs. (5,8) which gives opposite side correlations) estimating its 
contribution to two-photon correlation function in central rapidity region. 
Suggesting fast transition we can evaluate the contribution in the vicinity of 
fixed proper time r c . So one only has to shift to the rest frame of each element 
of the expanding volume and perform reintegration. Then the extention of 
the correlator in Eq.(7) to the case of expanding volume takes the form: 

2fc(a(ki)a(k 2 )) = G(k 1T + 

k 2T)4r 2) (21) 

with 

l|? = / £x T J dv ^^ rikx coshr?) (22) 

where we neglected n(k) in comparison with unity (see above). Therefore 
the normalized two-photon correlation function is given by (cf Eqs. (8-9)) 

C(k 1T , k 2T )\ yi=y2=0 = 1 + C HBT + R 2 (k T )G 2 (k 1T + k 2r ) (23) 
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with 



R(k T ) 




(24) 



Iqgp + 




We calculated the ratio R(kx) for different transition times 5r = Ofm/c, 
0.5/m/c, 1/m/c up to k T = 500MeV. In the region k T < WOMeV the ratio 
R is sizable for all these St being equal 4.94 at kr = 0, reaching maximal 
value R ~ 6 at A; T ~ 20MeV and falling down at k T = lOOMeV to i? = 1.78 
for 5t = 0, R = 0.95 for <5r = 0.5/m/c, i? = 0.55 for 5t = lfm/c . At 
larger transverse momenta the behaviour of the ratio R depends strongly on 
transition time 5r: in k T interval (200 — 500)^^6^ the ratio R(k T ) rises from 
1.3 to 3.9 for 5t = 0, it is approximately constant (R=0.3) for 5t = 0.5/m/c 
and it decreases from 0.15 to 0.03 for 5t = 1/m/c. On the whole it seems 
that the measurement of R(kx) gives a possibility to identify the transition 
effect, especially if 5r < 0.5/m/c 

5 Conclusion 

Estimation of photon emission accompaning transition between quark-gluon 
and hadron states of matter in heavy ion collisions shows that opposite side 
photon correlations can serve as a sign of the transition if transition time is 
small enough. 
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